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Abstract 



We propose a new dark matter candidate, "quirky dark matter," that is a scalar 
baryonic bound state of a new non-Abehan force that becomes strong below the elec- 
troweak scale. The bound state is made of chiral quirks: new fermions that transform 
under both the new strong force as well as in a chiral representation of the electroweak 
group, acquiring mass from the Higgs mechanism. Electric charge neutrality of the 
lightest baryon requires approximately degenerate quirk masses which also causes the 
charge radius of the bound state to be negligible. The abundance is determined by an 
asymmetry that is linked to the baryon and lepton numbers of the universe through 
electroweak sphalerons. Dark matter elastic scattering with nuclei proceeds through 
Higgs exchange as well as an electromagnetic polarizability operator which is just now 
being tested in direct detection experiments. A novel method to search for quirky dark 
matter is to look for a gamma-ray "dark line" spectroscopic feature in galaxy clusters 
that result from the quirky Lyman-alpha or quirky hyperfine transitions. Colliders 
are expected to dominantly produce quirky mesons, not quirky baryons, consequently 
large missing energy is not the primary collider signal of the physics associated with 
quirky dark matter. 



1 Introduction 



Dark matter is a big mystery. Weakly interacting massive particles (WIMPs) provide one 
interesting class of particles to serve as dark matter. There are, nevertheless, two main 
puzzles with typical WIMP candidates: 



(1) Abundance is determined by thermal freezeout. While a thermal freezeout origin can 
yield the observed cosmological abundance when the annihilation cross section is tuned 
to roughly 1 pb, this mechanism is entirely unrelated to the origin of matter, which 
arises from an asymmetry. The observational relation between the dark matter density 
and the baryonic density, po — 5pB, is a coincidence. 

(2) Elementary WIMPs with electroweak interaction strength couplings to standard model 
(SM) fermions generically have very strong constraints from direct detection bounds. 
To be "safe", WIMP interactions with the SM must be sub-weak strength, and typi- 
cally, their mass determined by a mechanism unrelated to electroweak symmetry break- 
ing. 



We propose a new model of dark matter that tackles both problems. The first problem can 
be addressed if the dark matter abundance is linked to the baryon abundance. This has been 
considered before, for example, in the context of technibaryon dark matter [1-6]. In these 
models, electroweak sphalerons re-process baryon and lepton asymmetries into technibaryon 
asymmetry. The constituents are electroweak charged, while the technibaryon composite 
dark matter is neutral. The sphalerons generate a relation between the number densities of 
leptons, baryons and technibaryons, 

ne - ni rih - ni ntb - n^i, (1.1) 

where n£ — ni, — n^, and rith — rifi, represent the lepton, baryon and technibaryon asym- 
metries, and the exact proportions are 0{1) depending on the electroweak charges of the 



technibaryon constituents. Cosmologically Pd/pb ~ 5, so that the relation Eq. (1.1) implies 
the dark matter mass M ~ 5 GeV. If this were the end of the story, technibaryon dark mat- 
ter (or any other weak scale model which connects the dark matter to the baryon asymmetry, 
e.g. [7,8]) would be ruled out. However, if the dark matter constituents are just becoming 
non-relativistic as the sphalerons are decoupling, there is an exponential Boltzmann suppres- 
sion in the technibaryon asymmetry relative to the lepton and baryon asymmetries, implying 
the dark matter can naturally have an electroweak scale mass. The other possible solution 
to the dark matter baryon coincidence places the GeV scale dark matter in a hidden sector 
weakly coupled to the SM sector [9-17]. 

The second problem can, in fact, be ingeniously solved by compositeness. In technicolor 
theories, technifermions interact with the SM through renormalizable interactions, while 
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a composite technibaryon dark matter candidate is charge- and electroweak-neutral. This 
ehminates renormahzable interactions with the SM below the electroweak scale, leaving only 
higher dimensional operators [3,4]. In the early 1990s this was thought to be an unfortunate 
result - dark matter could not be observed in the-then foreseeable future. Given the direct 
detection bounds today (e.g. [18,19]), this becomes a "feature" of our composite dark matter 
model that we exploit to naturally suppress the direct detection cross sections. 

The residual electroweak-mediated direct detection cross section of composite dark matter 
arises from "form factor" suppression, somewhat analogous to the suppression of the photon 
coupling to neutrons at energies much smaller than Aqcd- For example, the leading order 
operators that couple scalar technibaryon dark matter to the SM arise at dimcnsion-6 (charge 
radius) [4] and dimension-7 (chromomagnetic polarizability) [3], suppressed by two or three 
powers of Atc- 

What was not fully appreciated in the 1990s is that both of these operators can be elimi- 
nated. The charge radius vanishes in a limit in which the current masses of the constituents 
are equal. The chromomagnetic polarizability vanishes when the constituents do not carry 
ordinary QCD color. While this suggests rethinking technicolor dark matter (e.g. [5,6]), 
the model building difficulties of realizing a fully successful technicolor model incorporating 
flavor as well as avoiding electroweak precision constraints remains daunting. 

In this paper, we take a different approach, in the spirit of the Hidden Valley [20] and 
Kang and Luty's quirks [21]. We retain the new strong dynamics, but discard their role 
in electroweak symmetry breaking. The new strong dynamics gets strong at a scale below 
the electroweak scale. We call the candidate that arises in this approach Quirky Dark 
Matter (QDM). The simplest quirky dark matter model, as we will describe, contains exactly 
the same gauge and matter content as that of minimal SU(2) technicolor with two flavors. 
Amusingly, what was originally a problem of minimal technicolor - namely the possibility 
that the vacuum aligned to an electroweak-preserving state [22] - is now a "feature" here 
since we utilize the ordinary Higgs mechanism to break electroweak symmetry. Indeed, we 
do not want the strong dynamics to break electroweak symmetry (even a little bit) lest we 
run into electroweak precision constraints. 

With an ordinary Higgs present, quirks can obtain their mass through the Higgs mech- 
anism just like quarks and leptons. This has several implications: New contributions to 
the electroweak oblique corrections arise; we estimate them below. Quirky dark matter can 
interact with nuclei of direct detection experiments through Higgs exchange; this leads to 
an ordinary elastic scattering cross section that is right near the current bounds for a light 
Higgs boson. Finally, assuming the new strong force confines, new operators involving the 
Higgs are present that can allow the new glueballs to decay, independent of the quirk mass. 
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SU{2)q SU{2)l U{1)y U{1)qb 



2 2 +1/2 

2 - -1/2 -1/2 

2 - +1/2 -1/2 



Table 1: Particle content and charges under the gauge and global symmetries. 



2 Model 



2.1 Field Content 



The model of quirky dark matter that we mainly wish to consider consists of two flavors of 
quirks transforming under a new strongly interacting sector, SU{2)q, that hereafter we call 
"quirkcolor." Variations of this, with different numbers of quirk flavors and quirkcolors are 
also possible; we will remark on the possibility of more quirk flavors later in the paper. We 
assume that the quirkcolor coupling constant gets strong below the weak scale. The particle 
content and charges under the gauge and global symmetries are given in Table [l] in terms 
of two-component spinors. This assignment is chiral under the electroweak group, and thus 
requires Yukawa interactions with the Higgs, 



to give current masses to the quirks, rrig = XgV, for q = U,D. We enforce a global U{1)qb 
"quirky baryon number" that forbids the mass terms ^qC,q and ^.o^d and ensures our quirky 
dark matter candidate is stable (on at least cosmological timescales). 

Since QDM contains additional matter that acquires mass exclusively through the Higgs 
mechanism, there are new contributions to the electroweak precision parameters. The quirks 
in our model are weakly-coupled at the scale of their mass, and so we can perturbatively 
calculate AS, AT [23]. Generically, AS' = l/(37r) ~ 0.1, while AT depends on the splitting 
within the quirky electroweak doublets. As we will show below, eliminating the charge 
radius operator requires negligible splitting between the current masses of the quirks. As 
a consequence, the contribution to T from this sector is negligible. The minimal model 
therefore has AS" ~ 0.1, AT ~ 0.0, which is roughly at the 95% CL contour when comparing 
against LEP electroweak working group fits [24,25]. Suffice to say it is a straightforward 
(but unenlightening) exercise to slightly extend the model to give a additional contributions 
to T (and S) that result in a model fully consistent with electroweak precision data. 

In the minimal model, with the only particles transforming under SU{2)q given by that 
in Table [l| the quirkcolor group confines. The global S't/(4) ~ 5*0(6) symmetry is broken 
to S0{5) and we have 15 — 10 = 5 pseudo-Nambu-Goldstone bosons. The large current 
quirk masses ensure these composites are massive, forming "quirkonia" bound states with a 
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(2.1) 
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spectrum similar to heavy quarkonia [26]. We assume the quirkonia masses are sufficiently 
heavy to avoid LEP II bounds (i.e., quirks heavier than about 100 GeV). Existing Tevatron 
bounds, as well as prospects for future detection at Tevatron and LHC, will be studied in a 
future paper. 

Confinement of SU{2)q leads to quirkcolor glueballs. These glueballs decay through 
higher dimensional operators into SM matter. Depending on the scales, however, their 
lifetime may be very long [21,27], potentially leading to cosmological problems depending 
on the quirk masses and the confinement scale. In QDM, there are additional operators 
due to interactions with the Higgs. These interactions are written with an estimate of their 
contribution to the glueball width in Sec. [6j 

Our quirkcolor group SU{2)q, however, does not necessarily need to confine, if additional 
(massless) flavors transforming only under SU{2)q are present. This provides an interesting 
possibility in which quirkcolor flows to a conformal field theory without confinement. We 
emphasize that, for much of our discussion below, essentially none of our calculations depend 
on the scale (or existence) of confinement, so long as it is sufficiently smaller than the 
inverse Bohr radius of the bound states so that reliable non-relativistic calculations can 
be performed. To this end, we need the quirkcolor coupling, evaluated at the scale of the 
inverse Bohr radius, to be perturbative. The situation is quite analogous to heavy quarkonia. 
Indeed, we employ much of the formalism of non-relativistic effective theories developed for 
quarkonia and apply it directly to QDM. The systematic derivation of the non-relativistic 
limit from the relativistic degrees of freedom, following the quarkonia literature, is outlined 
in Appendix |X} 

Our composites include "mesons" and "baryons" depending on whether or not they 
carry the nonzero U{1)qb quirky baryon number. To satisfy the Pauli exclusion principle, 
the wavefunction of baryons must be anti-symmetrized with respect to all quantum numbers 
leading to a detailed spectrum of allowed states. Details outlining the construction of the 
bound state spectrum from the relativistic spinors to the non-relativistic mesons and baryons 
are given in Appendix |Xj Here we simply quote the results and present the baryon spectrum 
that is relevant for dark matter and its interactions, given in Figure [T} 

The dynamics and binding energies of the quirky baryons is our primary interest in this 
paper. For this, we need to construct the non-relativistic potential, V^(r). Formally, this 
is possible in the limit ^ rriqV ^ rnqv"^, where the potential is dominated by single 
quirkcolor gauge field exchange between the massive quirk^ Since our binding quirkcolor 
force is non-Abelian, the strength of the potential depends on the representation of the 
constituents. For our model, the quirks are in the fundamental representation, which gives 

^If SU{2)q confines at a scale Aq, we actually need mgV ^ m^u^ > Aq. This is because in the alternate 
Hmit TUgV » Aq » rnqv"^ one must integrate out the physical scale Aq before the potential can be properly 
defined. This procedure leads to an additional non-perturbative part in the potential that contains non-local 
quirkcolor gauge field correlators. 
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Figure 1: Sketch of the quantum mechanical energy spectrum of our quirky dark matter 
composite with the ground state and several excited states shown. Our notation Bg corre- 
sponds to baryonic states with total electric charge q and total spin S. We have included 
0{a^) (quirky Lyman-alpha) and 0{a^) (quirky hyperfine) splittings, but do not show the 
0(a^) (quirky Lamb shift) splittings or other higher-order effects. The lightest electrically 
charged baryons Bf (not shown), have spin one, and are slightly heavier than B^ due to 
subdominant electromagnetic corrections to the potential. 
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a potential 



Vs{r) = (2.2) 
r 

wher^ a(r) = 6*2(2) aQ(r) = |aQ(r), and Q;Q(r) is the strength of the quirkcolor force 
evaluated at the scale 1/r. This potential strictly arises from the one quirkcolor gluon 
exchange in the singlet channel, (^a^ij^ki)^ji — C2{2)6ij. This potential is similar to the 
Coulombic potential used to determine the bound states of the hydrogen atom. However, the 
non-Abelian nature of the quirkcolor binding force implies aqlr) is scale-dependent, which 
to leading-log is given by the /?-f unction [28,29] 

«Q(r) = aQirs) (l + ^^(H - Nf) ln(r/r^)) . (2.3) 

Here Tb = [a{rB)fj]~^ is the analogue of the Bohr radius in the hydrogen atom and gives 
the typical size of the bound states. We have written the r-dependent correction to the 
potential for a general number of quirkcolor flavors for completeness. The model on which 
we concentrate our attention has Nf = 2, as defined before in Table [T| 



2.2 Binding energies and splittings 

Our quirkcolor singlet bound states can be described by a Schrodinger-like equation with a 



quirkcolor force potential given by Eqs. (2.2), (2. 3). In the limit that the logr piece can be 



neglected, the potential becomes purely Coulombic - the same as the hydrogen atom - with 
the replacements mg ^ fi and ^ ot- One can formally approach the Coulombic limit 
if enough additional flavors are present to lead to a nearly scale-invariant quirkcolor theory 
while a remains perturbative. 

We will be interested in the regime where a is perturbative but not necessarily small, 
and with exactly two flavors of quirks as given in Table [Tj Hence, our non-relativistic 
potential has unavoidable logr dependence. In our calculations below, we express the effect 
of the log term as coefficients that multiply the exact solutions obtained in the case of a 
purely Coulombic potential. The coefficients have been computed by numerically solving 
the Schrodinger equation including the log term for a few specific choices of a. 

The binding energies of the n-th excited state of the quirkcolor singlets is given by 

En = (2.4) 

expressed in terms of the reduced mass of the quirks, = l/mu + l/mci. The n-dependent 
constant kn captures the difference between a pure Coulombic potential and our non-Abelian 



^We drop the subscript Q when writing a to keep the notation as simple as possible. 
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quirkcolor theory. Using our numerical evaluation of the Schrodinger equation, we find the 
energy levels of the first two states are corrected by ki ~ (1.2, 1.3, 1.4) and k2 — (1.9, 2.3, 2.8) 
for a{rB) = (0.2,0.3,0.4) and Nj = 2. 

Next, the hyperfine structure. The spin of each of the constituent quirks generates a 
dipole quirkcolor "magnetic" field which leads to a spin-spin interaction, 

Hhj = ^^^S,-S26%r) + ... , (2.5) 
3 murriD 

where the terms we have neglected do not contribute to the angular momentum £ = states. 
This is the leading non-relativistic contribution to the hyperfine structure. (Relativistic cor- 
rections, such as Thomas precession [30], are small so long as the quirk masses are much 
larger than the strong scale.) Sandwiching this Hamiltonian between states of the unper- 
turbed potential gives a splitting proportional to the (unperturbed) wavefunction at the 
origin, 

|V'i,o(0)| = Ci,o . (2.6) 

TT 

This is the familiar result from the hydrogen atom, except for the constant Cn,^ which differs 
from one due to the logr term in our non-relativistic potential. Numerically calculating 
the coefficient for (n,£) = (1,0), we find ci,o ^ (0.5,0.4,0.3) for airs) = (0.2,0.3,0.4) and 
Nf = 2. The hyperfine splitting is thus estimated from Eq. (2.6) to be 



^ ^ ^ /^^ «^ f 2 spin-triplet ,^ 

^^'°3mumD 1 -6 spin-singlet. ^ ' ^ 

As long as the quirk masses satisfy \mu — niul < E^f, the electrically neutral spin-singlet 
baryon is lighter than the electrically charged q = (-1-1,0, —1) spin-triplet baryons Bf, 
in agreement with [1]. This requires our quirk current masses to be very nearly degenerate, 
mu — Hereafter, we use Bq to denote our quirky dark matter candidate in its ground 
state, n = 1. We illustrate the spectrum of the ground and first excited baryonic states in 
Fig. HI 



3 Quirky Dark Matter Abundance 

Stable quirks transforming under a chiral representation of the electroweak group have an 
abundance that is necessarily related to the baryon and lepton abundance through the elec- 
troweak phase transition. That such a relationship is inevitable was discovered in early 
work on the technibaryon abundance from technicolor theories [2]. There it was shown that 
baryons and technibaryons could have a common origin, since sphalerons intermix baryon, 
lepton, and technibaryon numbers. More interestingly, Ref. [2] showed that the large mass of 
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the technibaryons causes an additional Boltzmann suppression of their abundance, roughly 
scaling as exp[— M^,/T*] where is the mass of the technibaryon at the critical tempera- 
ture where sphalerons shut off. This allows for TeV mass technibaryons to nevertheless 
yield roughly the right dark matter abundance today (for recent calculations in technicolor 
theories, see e.g. [5]). There are three global flavor quantum numbers - baryon, lepton, and 
technibaryon number - while the sphaleron violates only one linear combination, leaving two 
anomaly-free invariants [7]. Hence, while baryon and technibaryon numbers are related to 
one another, one cannot determine technibaryon number uniquely from just baryon number. 
Instead, baryon, lepton, and technibaryon numbers are ultimately determined in terms of 
linear combinations of two invariants which can be taken to be S — L number [31] and another 
combination involving both baryon (or lepton) number and technibaryon number [7]. 

The abundance of quirky dark matter is determined by an analysis similar to that of 
technibaryon dark matter. The main difference between our calculation below and that of [2] 
is that quirkcolor is assumed to be weakly-coupled through the electroweak phase transition. 
Sphalerons therefore yield an (asymmetric) abundance of quirks instead of quirky baryons. 
To also exploit the Boltzmann suppression of quirky baryon number, quirks must acquire 
mass before sphalerons shut off, which can occur if the electroweak phase transition is not 
first order. The Boltzmann suppression for the abundance of quirks is therefore proportional 
to exp[— Agf (T^,)/r^,], where v{T^) is the electroweak vev at the critical temperature T^. The 
final ratio of quirky dark matter abundance to baryon abundance is determined by three 
quantities: the two primordial anomaly-free U{1) invariants and the ratio mq{T^) /T^ = 
Agf(T*)/T*. In principle, f(T*) and T,, can be calculated within our theory. This requires 
incorporating the effects of quirks on the phase transitioi]|^ 

Following the classic calculation of [33] , the divergence of the baryon, lepton, and quirky 
baryon currents can be constructed from 

^ ~ 647r2 ^Mi'^P'^ 

where only SU{2)\y effects on N electroweak doublets are considered. It is straightforward 
to determine that the sphaleron carries B = Ng/2, L = Ng/2 and D = Nd/2 charge where 
A'g = 3 is the number of SM generations and Nd = 1 is the number of electroweak doublets 
that carry quirky baryon charge. We normalize the quirks to carry 1/Nq = 1/2 quirky 
baryon charge, given Nq = 2 quirkcolors, precisely analogous to the 1/Nc = 1/3 baryon 
number normalization of quarks. 

This result implies sphalerons violate the global U{1) number B + L + j^D. The or- 
thogonal combinations that are preserved can be taken to be Ji = B — L and I2 = B — 
(or L — [7]. Using the standard techniques [31,34], we enforce the following relations 

■^Examples of theories with additional chiral fermions with large Yukawa couplings have been considered, 
e.g., [32]. There it was found that the electroweak phase transition was weakened (without superpartners), 
which is not inconsistent with our expectations, though we leave a more detailed analysis to future work. 
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among the chemical potentials: (i) electric neutrality, and (ii) the vanishing of the chemical 
potential of the Higgs boson. With these conditions, and taking A^^^ = 3 and Njj = 1, we 
obtain 

o ^ m{x)+Af{xf]h + [n + 2f{x)]h , . 

17 + 113/(x) + 13/(x)2 ^ ■ ' 

J ^ -[17 + 77/(x) + 9/(x) + [17 + 2f{x)]h , . 

17+113/(x) + 13/(x)2 ^ ' 

n _ j36 + 4/(x)]/i-[lll + 13/(x)]/2 

^ 51 + 339/(x) + 39/(x)2 ^^'^^ 

where 

= ^/ r^—, (3.5) 

27r^ Jo cosh^ I ^Jy^ + 

in terms of x = mq{T^) /T^.. The mass density ratio is, therefore, 

PR = = f( ) [36 + 4/(x)]/i-[lll + 13/(x)]/2 M 

PB Brup ^^''^6/(x)[18 + 2/(x)]/i + [51 + 6/(x)]J2mp ' ^'^ 



where M is the mass of Bq in its ground state. In Fig. [2] we show contours of the resulting 
quirky dark matter density Pd/pb within the parameter space of the two primordial invari- 
ants Ii and J2. Pure leptogenesis, which corresponds to Ji = —L and I2 = 0, immediately 
implies D = B/3, independent of f{x). This mechanism is not viable since the mass of 
quirky dark matter would have to be M ~ 5mpB/D ~ 15 GeV, which is ruled out by direct 
collider searches. Pure baryogenesis or some mixture of all three remains perfectly viable so 
long as the lepton number of our universe remains unknown. 



4 Prospects for Direct Detection 
4.1 Overview 

There are three basic ways that quirky dark matter could potentially be found in direct 
detection experiments: 

(i) Elastic scattering through Higgs exchange. 

(ii) Elastic scattering through higher dimensional operators. 

(iii) Inelastic scattering through an excited quirky baryonic state. 
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B -L 



B -L 



Figure 2: Contour plots of the density ratio PdIpb = (1,5,25) shown by dashed, sohd, 
dot-dashed (red, blue, green) lines. The axes are the invariants (Ji, I2) = {B — L,B — 3D) in 
arbitrary units; a mirror symmetric plot can be obtained taking (/i,/2) {—h, —Ii)- Plot 
on the left has M = 200 GeV, x = 0.25, and on the right M = 1000 GeV, x = 0.25. The 
plots demonstrate that a viable region exists with Pd/pb — 5, corresponding to a "bathtub 
ring" around a valley in (/i, I2) space. The bottom of the valley has Pd/pb — 0. 

The third way, inelastic scattering, has been considered before in general [35] and recently 
in the context of composite inelastic dark matter [36,37]. Quirky dark matter is more akin 
to Ref. [35], where it was shown that one needs fairly small splittings, up to about 10 MeV, 
to allow for inelastic recombination. The smallest splitting in quirky dark matter is the 
hyperfine splitting. Combining a rough bound from LEP II, rriq > 100 GeV, with a ^0.1 
to satisfy direct detection bounds (explicitly shown later in this section), we find the hyperfine 
splitting Efif ~ 2a^p/3 > 30 MeV. So, we do not anticipate inelastic scattering or inelastic 
recombination in direct detection experiments. 



4.2 Higgs Exchange 

Our quirks acquire mass through the Higgs mechanism, and hence Bq has interactions with 
matter through Higgs exchange. Just as the Higgs couples to the qq content of the nucleon 
through {N\qq\N), the Higgs also couples to the quirk-quirk content of our quirky baryonic 
dark matter (B^lqqlBQ) . Unlike the nucleon, however, the quirkcolor gluon condensate is 
presumed to give a negligible contribution to the quirky baryon mass. The calculation of 
Higgs exchange is most easily done in the low energy effective theory below the scale of the 
quirky baryon. Then we can treat Bq as simply a complex scalar with a renormalizable 
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interaction with the Higgs, 

£ ~ MhB^o*B^o (4.1) 

where this interaction assumes the mass of Bq arises mostly from the current quirk masses, 
i.e., M ~ rriu + mjy. With this interaction, we can use the results of Ref. [38,39], which 
considered the scattering of real scalars through Higgs exchange, and read off the direct 
detection cross sectiorEl We obtain 

a(nucleon) = ^i^^(Z/, + (A - Z)/„)^ (4.2) 

where fj,{D,n) is the reduced mass of the Bq and nucleon for scattering off a nucleus with 
atomic number Z and nucleon number A. The nucleon parameters can be written as 



/. '^nucleoli 

/nucleon 



E/. (nucleon) 6 „ (nucleon) 

_q=u,d,s 

We have factored out the Higgs coupling to Bq, so that only nuclear physics-dependent 
coefficients are present. Numerically [40], 

/jfj = 0.020 = 0.026 (4.4) 

/^^ = 0.014 /g = 0.036 (4.5) 

and [41] 

q=u,d,s 

The strange quark contribution to the nucleon is much more uncertain. A recent lattice 
calculation suggests it is much smaller than has been previously estimated [42] (see also [43]), 
from which we extract f^g""^ — 0.02. 

Interestingly, since fi{D,n) ~ ^nucleon, there is essentially no dependence of the nucleon 
cross section on the mass of the dark matter. This occurs because the (mass)^ cancels 
out between the numerator (its coupling to the Higgs squared) and denominator (from the 
non-relativistic expansion of the cross section). Putting all of this together, we obtain 

a(nucleon) ~ 1.8 x 10"^=^ ( 111^^ cm^ . (4.7) 

\ fiih J 

The current best bounds come from CDMS [19], o"(nucleon) < 0.8-3.5 x 10~^^ cm^, for dark 
matter mass between about 200-1000 GeV. This means that if the Higgs is very near the 
LEP bound, quirky dark matter should be seen in direct detection experiments in the very 
near future. On the other hand, the absence of a direct detection signal would put a lower 
bound on the Higgs mass that is consistent with quirky dark matter. 

''^Note that their hS^ has a coefficient of Aw246 which translates into a Feynman rule with coefficient 
2Aw246 to account for identical particles. For a complex scalar, the equivalent Feynman rule constructed 



(4.3) 



from Eq.(4.1 1 has no factor of 2, and so our o'(nucleon) is smaller by a factor of 1/4. 
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4.3 Higher Dimensional Operators 



The interaction of quirky dark matter with nuclei in direct detection experiments can also 
proceed through higher dimensional operators involving the photon. Since Bq is an electri- 
cally neutral scalar composite, all its moments vanish. The leading interactions are then the 
charge radius and the polarizability operators [44], 

where is the charge radius, and Opoi is the electromagnetic polarizability of B^. These 
interactions do not scale with the mass {A number) of the nucleus, and so we cannot use the 
usual (T(nucleon) cross section to compare with experimental results. Instead, we derive the 
nuclear cross sections that result from the charge radius and polarizability. These nuclear 
cross sections are in principle easy to compare with experiments, except that experiments 
often quote bounds only on the average nudeon cross section. Assuming the mass of a detec- 
tor is dominated by one (heavy) isotope of a nucleus with atomic number A, the relationship 
between the nucleon and the nucleus elastic scattering cross sections are related by 

cT(Nucleus) = ^^7?f— cT(nucleon) , (4.9) 
IJ,[D,ny 

where /x(-D, n) and /x(-D, A^) are the reduced mass of the dark matter with the nucleon and 
the nucleus respectively. 



4.4 Charge Radius 

The leading order interaction of a photon to a neutral scalar bound state of charged con- 
stituents is given by the charge radius. The scattering cross section off a nucleus due to its 
charge radius is given by [44] 

(7(NucleUs)charge radius = — fx{D , Nf al^Z\% (4.10) 

To gain a feeling for the size of the existing constraint, for example from CDMS, we can 
compute the bound on r% for a few specific choices of dark matter mass: 

r (510GeV)-2 M = 200 GeV 
rl < } (440 GeV)-2 M = 400 GeV (4.11) 
[ (370 GeV) -2 M = 800 GeV. 

We now compute the charge radius of Bq in terms of the mass difference Sm^ = {mjj — mr,) 
and average mass rUq = {mu + mD)/2 of the quirks. The magnitude of the charge radius is 
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estimated by assuming the charge distribution inside the bound state takes the form 



p(r) 



{pu{r) - pDir)) 



qe 



V2 (mjtt)^''^ exp {—rriiar) 



(4.12) 



where pi = |(i?Q|'?/'j)| is the probabihty density of finding the i-th quirk in the bound state 
Bq. We approximated individual quirk wavefunctions by assuming the other quirk to be 
fixed in space and normahzed it such a way that the total probability of finding the quirk 
is 1/2. The charge radius, interpreted as a measure of the size of the bound state when 
probed electromagnetically, is defined classically as the second moment of the spatial charge 
distribution, 



D 



3q SrUq 



+ 



(4.13) 



where we have assumed 6mq <^ niq. Using the constraint from Eq. (4.11), we find 



rrin 



< 



25 GeV 



X 



6.4 X 10^3 
8.6 X 10^3 
1.2 X 10^2 



M 
M 
M 



200 GeV 
400 GeV 
800 GeV 



(4.14) 



Hence, the quirk masses must be very nearly degenerate to avoid generating an electro- 
magnetic charge radius that exceeds the direct detection bounds. Interestingly, we already 
required 6mq < Ehj, to ensure the lightest quirky baryon is electrically-neutral. We see that 
a self-consistent picture has emerged in which nearly or exactly degenerate quirks ensures 
both that the lightest quirky baryon is electrically neutral as well as a negligible electro- 
magnetic charge radius. In Appendix |B| we demonstrate that the vanishing of the charge 
radius can result from imposing an exact discrete symmetry, "UD-parity" , which enforces 
mu = mo- 



4.5 Polarizability 

Having discussed and discarded the charge radius operator, we now move on discuss the 
electromagnetic polarizability. The scattering cross section due to the polarizability operator 
is given by [44] 

144 oP' , 

a(Nucleus)poi ^ —p{D,NfZ^al^,^, (4.15) 

where the nucleus is assumed to be a sphere of homogeneous charge with radius tq = \fA x 
1.2 fm. To again gain a feeling for the constraint, for example from CDMS, we compute the 
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bound on ctpoi for a few specific choices of dark matter mass, 



[106 GeV)^^ M = 200 GeV 
«poi < { (95 GeV)"^ M = 400 GeV (4.16) 
(85 GeV)"^ M = 800 GeV . 



We now calculate the polarizability of -Bq. To proceed, we first quickly review the stan- 
dard polarizability calculation. If an electric field S is applied to the bound state in the z 
direction, it causes a perturbation to the Hamiltonian 



pert 



qeSz 



(4.17) 



where qe is the constituent quirk charge. The Schrodinger equation in the presence of this 
perturbation. 



can be approximately solved by standard perturbation theory methods, 

{km 



\0)+qeSj2 



k>0 



En — El 



■\k) 



(4.18) 



(4.19) 



where Hoik) = Ek\k), and the ground state is |0). In a moment we will identify the states 
with the relevant states of Bq. At leading order, the dipole moment of this state in 
the z direction is composed of two terms: 



Pz 



-qe{'ip\z\ 



-qe\mQ) + Y, 



k>0 >- 



^^'^^^"'°^(0|.|A:)+ yf\ fc|.|0) 



En — Eh 



En — El, 



(4.20) 



The first term, (0|z|0), is the moment of the unperturbed state (if any). The second set of 
terms represent the moment induced by the electric field, Pind = ctpoiS, where Opoi is defined 
as the polarizability. We thus find 



O^pol 



2gV5:^ 



fe>0 



Eq — Ek 



(4.21) 



In order to extract the proper dependence of the matrix element in Eq. (4.21) on aem 
and a, we must also include the ordinary electromagnetic Coulomb potential. The un- 
perturbed potential is then simply V{r) = — {a{r) + q^aem) A, the effective Bohr radius 
then is given by = [{a{rB) + q^aem)lA~^ ^^"^ energy eigenvalues are given as Ej^ = 
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— {a{rB) + q^cxemT ■ Denoting the states by the usual quantum numbers, {\k)} 
{|n, /,m)}, with the ground state |0) = |1,0,0), we find 



a 



pol 



\{n, l,0|:|l,0,0)f 



n>l ^ " 

Koi I ^ r| (4.22) 



2 a + q^a, 



em 



The result above shows that there is an additional q^a^m/ + q^cxem) suppression in the 
polarizability relative to the hydrogen atom. This agrees with the analogous calculation for 
the electromagnetic polarizability of heavy quarkonia, substituting the quirkcolor coupling 
with the QCD coupling [45]. This factor arises since the binding potential is proportional 
to a + q^Uem- The additional non-Abelian correction resulting from the logr term in the 
potential is encoded in the coefficient kpoi- Numerically solving the Schrodinger equation, 
we calculated the corrections to the first few terms, finding the largest correction to the 
(2, 1, 0|2;|1, 0, 0) term and a smaller correction to (3, 1, 0|2;|1, 0, 0) term. Extrapolating from 
these results, our numerical estimate is kpoi = (1.3,1.4,1.5) for a^rs) = (0.2,0.3,0.4) and 
Nf = 2. 

We can now easily determine the bounds and prospects for detection of quirky dark 



matter as a function of a using our expression for apoi in Eq. (4.22 ), In Fig. |3jwe show theory 
predictions and experimental bounds on an effective cross section. To combine results from 
several experiments with theory predictions the effective nuclear cross section applicable to 
the polarizability operator is 

cTpoi,eff = ^ ^^^ a(Nucleus)poi . (4.23) 

where the individual experiments' target size (tq) and charge {Z value) have been factored 
out, along with two powers of the reduced mass, analogous to Eq. ( 4.9[ ). 



To gain an appreciation of the constraint on r^^, we can compute the lower bounds on the 



inverse Bohr radius derived from the polarizability constraint directly from Eqs. (4.22 ), (4.16 ): 




- 32 GeV M = 200 GeV 
r^^ > \ 33 - 27 GeV M = 400 GeV (4.24) 

- 24 GeV M = 800 GeV , 

where the range in value corresponds to 0.2 < a < 0.6. Since the cross section induced by the 
electromagnetic polarizability is proportional to r^, there is a relatively minor dependence 
on the mass of dark matter and the strength of the strong force. 
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Figure 3: Bounds of nuclear cross section. We have plotted a"poi,eff from Eq. (4.23) us- 
ing data from CDMS (blue dash), Xenon (green dot-dash) and Tungsten (red dot). The 
solid lines are theory predictions due to the polarizability operator for two flavors as a 
function of dark matter mass M/j. From right to left we have plotted crgg for a{rB) = 
{0.2 (blue), 0.4(purple), 0.6 (yellow)} respectively. 



5 Indirect Detection through Absorption Lines 



An interesting feature of quirky dark matter is that it has a spectrum of excited bound states 
that can be excited by absorption of SM particles such as the photon and Z. The photon 
interaction, in particular, provides an fascinating possibility to probe dark matter directly 
through photon absorption lines, entirely analogous to how matter itself is probed through 
its own photon absorption lines. The typical energies for quirky excitations, as we see below, 
are in the gamma-ray region. The possible existence of dark lines due to inelastic transitions 
of dark matter was considered before in a somewhat different setup in Ref. [46]. We will 
apply their results, suitably modified for our case (with one correction to their formulae 
given in Appendix [C]), to the two transitions of greatest interest: the quirky Lyman-alpha 
transition, and the quirky hyperfine transition, illustrated in Fig. [T] 
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5.1 Quirky Lyman- alpha 



The photon energy corresponding to quirky Lyman-alpha transition, between the ground 
state and the first excited state can be read off from Eq. (2.4), 

ELa = E2-Ei= ' -aV- (5.1) 

The deviation from purely Coulombic is characterized by the constant {Aki — k2)/3 ~ 
(1.0, 1.0, 0.9) for a{rB) = (0.2, 0.3, 0.4), which is neghgible for the quirky Lyman-alpha tran- 
sition. 

The width of this absorption feature can be computed in the Coulombic approximation 
by analogy with hydrogen Lyman-alpha. We obtain 

4 64 

TLa = -q^aemEU{l\f\0)\' = ^«em«V (5.2) 
9 6561 

This width, taking a —>■ and q —>■ 1 reproduces the well known value for hydrogen 
Lyman-alpha. Plugging in representative values, the width is 

^ 5xlO-Gev(Ay'(^) . (5.3) 

It will also prove convenient to express the width as 

-aemtt^ ^ 1.9 X 10"^ (5.4) 



^ La 512 2 -I ^ 1 n— 4 -2 



Eia 19683 



5.2 Quirky Hyperfine 

The next transition of interest arises from the quirky hyperfine structure. The energy split- 
ting between the and Bq states is given by 

2 

Ehf = Ci,o-aV (5.5) 
where ci^ is an order one constant that characterizes the deviation from purely Coulombic. 



A few specific values were computed in Sec. 2.2 



The decay rate of the quirky hyperfine transition can also be computed. The result in 
the purely Coulombic approximation is the same as for quarkonia [47] 



4 .E^ 



3 



aemq = ^OiemOi /i . (5.6) 

yU"^ 81 
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Plugging in representative values, we find the width to be 



r„ ^ 8xlO-Gev(A)"(_Ji_) (5.7) 



It will again prove convenient to express the width as 



^ = ^aema^ ^ l.lxlO-^a« (5.8) 



5.3 Absorption and Broadening 

For these transitions to be visible to gamma-ray observatories, three things are required: 
(1) the cross section on resonance should be substantial, (2) the width of the dark matter 
Doppler-broadened resonance should be resolvable, and (3) the energy of the transition 
should be within the range that is observable. Several observatories view the universe in 
the relevant energy range, satisfying (3). In this section we compute the cross section and 
Doppler-broadened width. 

The absorption cross section near a resonance takes the usual Breit-Wigner form [48] 

where M is the mass of the ground state (the mass of quirky dark matter), and = M + E^ 
is the mass of the excited state. The dominant decay of the resonance is into photons, 
r* ~ TMt->M-ry since decay to an on-shell Z is kinematically forbidden (throughout our 
parameter space), while decay through a virtual photon to a fermion pair is suppressed by 
aem- The factor of comes from 4^(2 jji + 1)/ ((2si + 1) (2s2 + 1) where Jr — lis the massive 
spin-1 resonance, (2si + 1) = 2 for the massless photon, and S2 = for the massive scalar 
quirky dark matter particle. 

There are two potential methods to exploit this absorption cross section. The first, and 
most promising, is to look towards massive galaxy clusters that have the largest column 
density of dark matter and a large dark matter velocity dispersion. We then compare this to 
the seemingly less promising method of looking for absorption over cosmological distances. 

Wc follow the formalism of [46] to determine the effect of the broadening and the size of 
the absorption cross section. This formalism applies to any of the photon-induced baryonic 
excitations, including quirky Lyman-alpha and quirky hyperfine. The resonance energy is 

= (5.10) 
which is roughly equal to the mass difference, — M for a ^ 1. 
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The dark matter velocity distribution within a galaxy cluster is assumed to follow a 
Maxwell-Boltzmann velocity distribution, denoted /mb(p)- The effective absorption cross 
section for a photon of energy E-^ is 



a{E,) = / rfp/MB(p)(a), (5.11) 
Jo 

where 

^ f' d cose Qn M^Tl 

^ ^ 7-12 pL(^-M2)2 + M2rr ^ • ^ 

is the total cross section after integrating over the incident angle. The center of mass energy 
is s = + 2Ej ^a/p^ + — j9cos^^ while the momentum in the center-of-mass frame is 
given by 

pL = ■ (5-13) 



The integral in Eq. (5.12) was solved analytically in the Appendix of Ref. [46], which we 
checked and agree with except for a correction to one expression given in our Appendix [C] 

For our case, with ~ M ^ E^, the analytic formula can be drastically simplified in 
the following limits. First, observe that M^kF* <^ AM^ which is equivalent to 2F* <^ E'!^'^^, is 
automatic for our perturbative calculation of the width of the resonance to be valid. Next, 
consider the limit 

F,M, < 2E;> (5.14) 

which corresponds to roughly 

F, < 2^^^V, (5.15) 

given the average momentum of dark matter with a Maxwell-Boltzmann distribution is 
roughly p ~ Ma^. Now compare this expression with Eqs. (5. 4), (5. 8). For a < 1, Eq. (5.15) 
is satisfied for > 10""^, which is itself satisfied by large galaxy clusters. 

Putting all this together, we obtain the following simple formula for the cross section on 
resonance, 

which can be integrated over a Maxwell-Boltzmann distribution to become 

aiE^ = = ^=7^3/2^^ . (5.17) 



'V 
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Interestingly, this form of the on-resonance cross section has several important features: 
First, there is no explicit dependence on the mass of the particle. Second, the resonance 
photon energy dependence in (from Lyman-alpha or hyperfine) exactly cancels the de- 
pendence in the denominator. This leads to very simple expressions for the on-resonance 
absorption cross section 

= 16384^2^ 4 
^ ^ ^ ' 59049 ^""(7^ ^ ' ^ ' 

7^3/2 2 

a{E,^^^^) = ae„.a'-^ {hf) (5.19) 

The opacity to 7-rays due to these dark lines on resonance can be estimated using the 
optical depth r = (7Edm/-/W^ where Edm is the dark matter surface density associated with 
the integral along the line of sight of the dark matter density. As an example, consider the 
Coma cluster, which has a surface density that was estimated by Ref. [46] to be S^m — 
5 X 10^^ GeV/cm^ with velocity dispersion = 820 km/s. Plugging in these characteristic 
values, 

, ^ ^ .820 km/s /25GeVy 200 GeV Sdm , 

T res = 1 X 10 ° — ; — — ^ (La) (5.20) 

' a, V / M 5 X 1026 GeV/cm' ^ ^ v ; 

, ^ ._2 820 km/s /25GeVV 200 GeV Sdm r^.J^^rr^.^ 

rlres = 1.5 X 10 M in26n../ 2 (^/) (5-21) 

(^v \ fB / iW 5 X 1026 GeV/cm 

we see that r is much smaller than one. For the Coma cluster, only a small fraction of 
photons are expected to be absorbed in either the Lyman-alpha or hyperfine transitions. 

There are three potential ways to improve on this result. The first is to search for systems 
with larger surface mass densities. This is most likely to arise in the largest clusters that are 
also the most compact, and thus have a very concentrated mass function. The second is to 
perform more precise measurements of the photon flux of particular clusters, which would 
allow probing optical depths considerably smaller than one. Third, combining 7-ray spectra 
from many different clusters of different mass, velocity dispersion, and redshift, and suitably 
processing them into a normalizable spectrum, one could significantly improve the search 
for dark lines through large scale galaxy cluster sampling. Given that the cluster number 
density is rising rapidly as the cluster mass is decreased (see, e.g., [49]), and with improved 
photon flux and energy resolution, this is probably the best approach for the future. 

Assuming methods are developed to effectively probe these small optical depths, it is 
also important to know the Doppler-broadened resonance width. In our case, the Doppler- 
broadened width arises from the dark matter velocity distribution in the observed systems. 
For a Maxwell-Boltzmann velocity distribution, we find the observed width of the resonance 
is well-fit to an Gaussian 



a{E^) = a(E^^^)exp 



(5.22) 
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with a fractional width lS.EjE ~ a^. 

The velocity dispersion of dark matter in galaxy clusters has been found to scale as [50] 

..=.(1080 km/=)(^^|^)"" (5.23) 

where h{z) = H{z) /lOO kms~^ Mpc~^ is the normalized Hubble parameter at redshift z and 
Mcius is the cluster mass (defined as the mass within a sphere encompassing a mean mass 
density of 200 times the background matter density at that redshift). Searching for dark 
lines of dark matter in clusters thus requires the gamma-ray energy resolution comparable 
to the velocity dispersion. The full width at half maximum (FWHM) is approximately 
2.35a„ ~ 0.003 — 0.01 for the largest galaxy clusters with mass between about 10^^~^^Mq. 
Interestingly, for soft 7-rays up to 8 MeV, this resolution was achieved with the INTEGRAL 
spectrometer [51]. The EGRET and Fermi/GLAST observatories extend up to much higher 
gamma-ray energies, 30 GeV and about 1 TeV, respectively. Unfortunately, the FWHM 
energy resolution of these instruments is roughly 0.2 for EGRET and between 0.05 — 0.1 for 
Fermi/GLAST, which is a just bit too course to likely resolve the 7-ray dark line feature 
that we predict in our model. 

Finally, we consider the effect of scattering over cosmological distances. At a redshift of 
2; ~ 1 or larger, a photon with a given initial energy will sweep out a resonance of width 
F ~ zE^ (e.g. [52,53]). In order for this to give an observable effect, however, the mean 
free path i for photon absorption must be shorter than the cosmological distance that the 
photon spends on resonance 

' - ^ < H-'^ . (5.24) 



Pd{o')\tcs ~ 

Dark matter in the cosmos has negligible momentum, and so the resonance cross section can 



be obtained from Eq. (5.9) in the limit p <^ M, (o")|res = Q'^/plm = Svr/i?^. Substituting, we 
obtain 

P , M f E^ Y , , 



H-^ 200 GeV V 100 MeV, 

which shows that even if F ~ i^..^, the photon does not travel nearly far enough to be absorbed 
over a cosmological distance. 



6 Quirkcolor Glueball Decay 

With exactly two flavors, our quirkcolor theory confines. It is straightforward to estimate the 
confinement scale Ag. Quirky theories suffer from a potential cosmological problem, namely 
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photon injection during nucleosynthesis, if glueballs decay into photons with a hfetime that 
is of order 1 s. 



Confinement gives mass to the quirkcolor "glueballs" that we crudely approximate to 
have mass Aq. For theories with vector-like quirks, Refs. [21,27] showed that the glueballs 
decay slowly, since the leading operators are suppressed by many powers of the quirk mass. 
One such operator is 

where Fq and F are the field-strengths of the quirkcolor group and electromagnetism, re- 
spectively. This results in a decay rate 

^ ~ ^ 8^^ ^g " [iGeV ' ■ ^^-^^ 

quirks ^ / 

where the sum is over all of the quirks given in Table [T} 

A different operator exists in our model due to the Higgs coupling to our chiral quirks. 
Integrating out quirks and the scalar Higgs boson simultaneously results in a dimension-7 
operator 

-^i^^Q ^Q.J f (6.3) 

where v = 174 GeV. This leads to glueball decay into a pair of light SM fermions that 
satisfies 2m j < Aq. Despite the lower dimensionality, this operator is not obviously less 



suppressed than Eq. (6.1), due to the Yukawa suppression ruf/v. The decay rate is 



^ ^ a{mg)^m} A^ ^ / Aq V / 115 GeV V fa^V / rrif y 
^ Stt v^ml VlGeVy V J \ O.l J VO.l GeVy ^ 



(6.4) 

We see this decay rate is comparable to the rate into photons for the example parameters. 
The main distinction we draw is that the Higgs-mediated decay does not depend on the quirk 
mass. Hence, we can contemplate quirk masses that significantly exceed 100 GeV without 
necessarily leading to cosmological difficulties of late decaying glueballs, so long as the Higgs 
is relatively light. We emphasize that while our estimates are parametrically correct, they 
nevertheless have significant uncertainties, particularly with regard to the identification of 
the glueball mass with Aq. 



7 Discussion 



We have presented a new theory of dark matter that is made up of a baryonic composite 
of a quirks that transform under a new strongly coupled sector, SU{2)q quirkcolor. The 
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abundance of quirky dark matter is linked to the baryonic abundance through electroweak 
sphalerons. 

The baryonic excitation spectrum was computed, including the fine and hyperfine struc- 
ture. The lightest baryonic state can be made automatically charged-neutral when the quirks 
have (nearly) degenerate masses, which can be ensured through a discrete symmetry (UD- 
parity). Degenerate quirk masses also eliminates the dimension-6 electromagnetic charge 
radius operator, allowing a much larger range of scales to be considered. 

Quirky dark matter is at the threshold of direct direction though elastic nuclear recoil. 
Two processes lead to nuclear recoil cross sections through (i) Higgs exchange, which couples 
proportional to the (mass)^ of the nuclei (as usual), and (ii) electromagnetic polarizability, 
which couples to the electric charge Z^. 

Indirect detection may be possible by searching for a gamma-ray "dark line" spectroscopic 
feature in galaxy clusters that results from the quirky Lyman-alpha or quirky hyperfine 
transitions. This is a difficult measurement that might be possible in the future. It requires 
sensitivity to optical depths much smaller than one. We envision this could be accomplished 
with excellent gamma-ray spectral sensitivity applied to a large number of galaxy clusters, 
suitably combining the results together. The feature itself has a Doppler-broadened FWHM 
roughly of order 2.35AE^/ ~ 2.35(j^ ^ 0.003 — 0.01. This is close to but somewhat smaller 
than the FWHM energy resolution of EGRET and Fermi. Indirect detection through other 
means, such as annihilation in the Sun, galaxy, or beyond, does not occur so long as the full 
global U{1)qb quirky baryon number is exact. Annihilation signals would reappear if quirky 
baryon number were broken to a Z2, and signals of dark matter decay would result if [/(1)qb 
were completely broken (but only very, very slightly). 

The collider signals of quirky dark matter represent a plethora of possibilities [21]. Quirks 
can be pair-produced, which form bound states with quirkcolor strings attached. They will 
stretch and fiop back and forth shedding angular momentum in some combination of quirky 
glueballs (from the quirks and the quirkcolor string) as well as photon emission from the 
quirks which may result in interesting underlying event signals [54]. Eventually the quirks 
bound in a mesonic state recombine and annihilate, somewhat analogous to heavy quarkonia 
annihilation. Quirky baryon production is kinematically suppressed due to the need to pair 
produce a baryon and anti-baryon (four quirks total), to conserve quirky baryon number. 
Clearly, the collider physics of quirky dark matter is an area ripe for future exploration. 



Appendix A: Non-relativistic Details 



The construction of the non-relativistic theory starting from the theory described in Sec. 2.1 



begins with the chiral two component spinors in Table. [T} Below the electroweak scale, the 
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two- component fermions can be written in terms of 2 four-component Dirac fermions and 
their charge conjugates for the quirks and the anti-quirks respectively, 



i 

where i runs over the flavor indices {U,D}. In this basis, the four component Dirac spinors 
constructed are in the Weyl basis, which makes chirality manifest and is the most convenient 
choice for representing a relativistic chiral theory. In a non-relativistic theory, however, the 
Dirac basis is more suitable. The 7 matrices and all the four-component spinors can be 
rotated from the Weyl basis to the Dirac basis by the transformation 

7m Ujf^U^ and /, ^ Uf, , where U=^ (^}^ . (A.2) 

The advantage of this basis is that 70 = diag(l, 1, 1, 1) is diagonal, and so the dominant 
component of the Dirac four spinor (namely, (1 + 70) fi), is a two component Pauli spinor, 
which furnishes the minimal representation for the non-relativistic spinor field. 

The next step is to eliminate the quirk mass scales rrii while keeping the heavy quirk 
fields. This is accomplished by a simple time-dependent rescaling of the Dirac fermions in 



Eq. (A.2). There is now a preferred frame, namely the center of mass frame, which is the 
frame in which we work from now on. Hence, the dominant component of the full Dirac 
spinors become 

= e^"^'* ^ (1 + 7o) /. (A.3) 

X^ = e-^'"'* ^ (1 + 70) /r • (A.4) 

The field ip annihilates a heavy quirk field, while x creates a heavy anti-quirk field. These 
spinors roughly correspond to particle and antiparticle and are appropriate for a non- 
relativistic approximation about the center-of-mass frame of reference. Classically, the entire 



quirk and the anti-quirk Dirac spinors in Eq. (A.l) can be written in terms of the spinors in 



Eqs. (A.3) and (A.4): 




f^ = e-^'"^* I .,.3 I and = e^^^* I , J 1 (A.5) 



2ini+iD(j 



Xi 



The non-relativistic Lagrangian is written in terms of the spinors i}) and x that desig- 
nate almost on-shell quirks and anti-quirks. The Lagrangian can now be computed as an 
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expansion in l/rrij [55]: 



{d-E,-E,- d) + I^V'l'? ■ X - X (A.6) 



1 



where iDq = ido — gqAg^^, iD = iV + igqAq. The electric and magnetic quirky gauge fields 
are defined as usual Eq. = Fq^- and Bq. = eijkF^^. At tree level, the matching is simply 
cf = cd = Cs = 1. This is corrected due to quantum effects by Oi^g"^). The terms in this 
Lagrangian in Eq. (A.6) have well known physical interpretation. The a • Bq term is the 
quirkcolor-magnetic moment interaction, the D ■ Eq term is the Darwin term and the D x Eq 
term is the spin-orbit coupling. The term contains the first relativistic correction. 



A.l Spectrum 

The ground state spectrum of baryons and mesons can be determined by cataloging the 
Lorentz invariant bilinears made out of two ip^s or one ip and one respectively. The 
decomposition of the fermion bilinears in (fiavor,spin) space can be written as 

(2, 2)® (2, 2) = (1„1J©(3„1J©(1„,3,)©(3„3,) , (A.7) 

where the subscripts s and a denote the symmetric and antisymmetric representation, re- 
spectively. Mesons can be written in all of these representations since they are formed of 
non-identical particles. Baryons, however, must satisfy the Pauli exclusion principle. Com- 
bining two identical quirks into a Lorentz invariant requires an antisymmetric SU{2)q con- 
traction, and consequently, only totally symmetric combinations of (fiavor,spin) are possible, 
specifically, la) and (3^, 3^). The resulting baryonic states can be written as 

(1„1,) : 5° = (A.8) 

(3„3.): = ^'"-{cT^'iaX^,, (A.9) 

in terms of ijj^"' = e^^e'^^ipj^. Here, fiavor indices are designated by Latin letters etc.), 
spin by Greek indices (a, P, etc.), and quirkcolor indices have been suppressed. The resulting 
states consist of: Bq, a complex scalar with zero electric charge and +1 baryon charge, and 
B^^ , a massive spin-1 vector that is triplet under fiavor and also carries baryon charge. 
The anti-baryons with opposite baryon charge are similarly constructed. The B^^ can be 
decomposed in terms of spin-1 baryons with electric charge q, denoted by Bl (suppressing 
the spin index S), as 

= (A. 10) 

B^ = (A.ii) 
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The mesons are made of x ^"^^ 4' iii representations given by Eq. (A.7). In rough analogy 
to mesons in QCD with two flavors, we can write the mesons as 

(A.12) 
(A.13) 
(A. 14) 
(A.15) 

We have suppressed all fermion indices. The proper combinations of flavor states that yield 
definite electric charge states are formed analogously to Eq. (A. 11). 



la, la) : 


V 


= 


3s, la) • 


TT 




Is, 3a) : 





= i^(T^x 


^3s, 3^) : 


P 


= ^cj^cj^X 



A. 2 Non-relativistic bound states in Quantum Quirky Dynamics 

We assume that the quirkcolor force is weakly coupled and thus described by perturbative 
physics with the quirks much heavier than the resulting bound state energies. In this limit, 
the timescale of the relative heavy quirk movement becomes much larger than the timescale 
of the quirky gluon dynamics. Then, feedback effects of the moving heavy quirks on gluons 
can be neglected, and so the adiabatic approximation or the leading Born-Oppenheimer 
approximation should be applicable. Also, one can use the different energy scales just like in 
the positronium problem in QED (discarding the annihilation effect). The non-relativistic 
bound state is characterized by the scale of the quirk mass m (hard), the scale of the 
momentum transfer p ~ mv (soft) and the scale of kinetic energy of the quirks in the center 
of mass frame E ~ /m ~ mv"^ (ultrasoft), where v is the heavy quirk velocity in this 
frame. In our weakly coupled non-relativistic system, f ~ -C 1 and it follows that the 
three relevant scales are hierarchically ordered {i.e. m ^ mv ^ mv'^). 

Interestingly, such a situation also happens in quarkonium physics, where the same hier- 
archy has been utilized to construct equivalent effective theories to describe the quarkonium 
spectra and interactions (for some reviews, see e.g. [56-64]). First, the transition is from 
the EFT with relativistic quarks to the non-relativistic effective theory (NREFT) with the 
quarks (and larger momenta) integrated out. The theory describes dynamics of heavy quirk- 
antiquirk pairs at energy scales in the center-of-mass frame much smaller than their masses. 



In Eq. (A. 6), we have reproduced the NREFT Lagrangian. In quarkonium physics, a higher 
degree of simplification has been achieved by exploiting mv ^ mv"^ and building the so-called 
potential-NREFT (or pNREFT) [56,58], where degrees of freedom of ~ mv are integrated 
out. In this way, an analytical calculation of the spectrum becomes possible. The ultra-soft 
degrees of freedom that remain dynamical in this theory are quirks of momentum mv and 
energy mv"^ and quirky gauge fields of momentum and energy less than mv"^. The matching 
of the pNREFT to the NREFT is perturbative as long as mf ^ > Aq. 

The low energy theory is described in terms of quirk bilinears, which depend on the 
relative distance between the two quirks, r = xi — X2, and the center-of-mass coordinate 
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— * 

it! = {xi -\- X2)/2. All gauge fields are multipole expanded in r, and therefore depend only 
on R. At leading order in the multipole expansion [57, 58] , 

C = S^(^dt + ^dl-Vs{r)^S . (A.16) 

In the above, S is any quirkcolor singlet field and K(r) is the matching potential. The 
reduced mass of the system is denoted by /x. For our model, 

i = ^ + ^. (A.17) 

Hence, at leading order in the multipole expansion, the equation of motion of the singlet field 
is simply the Schrodinger equation! Determining the bound state energies is thus very similar 
to a standard quantum mechanical calculation. The main difference is that the potential 
depends on a scale-dependent quirkcolor coupling that introduces logr dependence in the 
potential. 



Appendix B: UD-Parity 

In this appendix we demonstrate the charge radius operator is odd under a certain Z2 symme- 
try under which the [/-quirk transforms to D-quirk, which we call "UD-parity" . Imposing 
that the ultraviolet theory is UD-parity symmetric therefore eliminates the charge radius 
operator and automatically ensures the lightest baryon is electrically neutral. 

Consider the limit m^/ = m/j. In the ultraviolet theory, one can show that the Lagrangian 
is symmetric under the following transformation: 

Cc/^Cd and ^ ^ -^d and Xu ^ Xd (vk ^^^\ 

A^^A^ and A^' ^ -A^" ^ ' 

In the above and are the quirkcolor gauge fields and the electromagnetic gauge fields 
respectively. Under UD-parity, the electrically neutral scalar baryon is odd and the vector 
baryon is even. 

Bo -Bo B^^ B^^ (B.19) 

and thus 

Ei^-Ei Bi^-Bi rl^-rl. (B.20) 
Hence, r|) must vanish if the low energy theory preserves UD-parity. 
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Appendix C: Integration Results 



We have verified the results in [46] including their Appendix, except for the expression for 
their Cg in (A. 4). We find the correct expression is 

where AM^ = M^—M"^. With this correction, we were able to reproduce Ref. [46] 's numerical 
results (including their Fig. 3). 
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